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We extend the nonlinear spin wave theory (NLSWT) for the spin 1/2 antiferromagnetic Heisenberg 
model on a triangular lattice (TAFHM). This novel NLSWT considers the corrections one order 
higher in 1/5 than the linear spin wave theory (LSWT). It also distinguishes in which circumstance 
the negative energy excitation, the sign of the breakdown of LSWT, shall be renormalized to be 
positive both by a boson normal ordering and a self-consistent iteration. We draw a phase diagram 
by testing the stability of various magnetic orders for different parameters. In particular, the 
incommensurate configuration is found unstable by our study. The new phase transition point 
(PTP) of the collinear configuration agrees well with various previous studies. 
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I. INTRODUCTION 

5=1/2 TAFHM is extensively studied. This model 
may describe the low-energy physics of various quasi- 
2D triangular materials such as Cs 2 CuCLfi, k-(BEDT- 
TTF) 2 Cu 2 (CN) 3 2 , and EtMe 3 Sb[Pd(dmit) 2 ]2^ for which 
experiments have indicated a magnetically disordered 
Mott insulator, i.e., a possible spin liquid state, at low 
temperatures. Various theories are applied to understand 
how the geometric frustration and low dimensionality af- 
fect the physics. The diagonalization method^ is exact 
but limited in the lattice size, while the finite size ef- 
fect is so strong to even change the phase diagram 5 . 
There are also nearly-exact numerical simulations such 
as DMRG-. They have encountered problems with over- 
whelmingly large entanglement when the dimensionality 
increases to 2D, although our own numerical work using 
the Entanglement Perturbation Theory (EPT), moves a 
small step further in this direction. The numerical sim- 
ulation is powerful as we already see very rich details 
of the physics which is believed to be closely related to 
the geometric frustration, low dimensionality and strong 
quantum fluctuation. An example is the phase change 
due to the finite-size effect mentioned above. 

However, currently, the most efficient theoretical 
framework to understand the model in the thermody- 
namic limit is the spin wave theory (SWT), which so far 
unfortunately only gives approximate description. There 
are two well known questions on the validity of SWT, i.e., 
if the assumed magnetic order is physical and if its ap- 
proximation is sufficient when the first answer is positive, 
because the conversion to the bosonic Hamiltonian in an 
expansion of Holstein-Primakoff transformation 7 restores 
the spin statistics only either in large S limit or when 
higher order corrections are properly taken into account. 
To get an insight into the answers, one usually examines 
the excitation spectrum. First, there must be Goldstone 
modes in the K-space to reflect the symmetry of both the 



lattice and the assumed order. It is automatically satis- 
fied in most SWT including LSWT because the unitary 
transformation, that is the starting point of SWT, natu- 
rally inherits those symmetries. Special cares, however, 
should be taken if it is involved with a self-consistent 
iteration or when there is a singularity cancelling issue 
in the perturbation. Second, the excitation spectrum 
must be positive. The emergence of the negative energy 
is a sign of the overwhelming quantum fluctuation over 
the assumed classical configuration, implying the break- 
down of SWT. Furthermore the interplay of the above 
two factors can sometimes lead to subtle situations. Par- 
ticularly for a spin 1/2 TAFHM, LSWT has a series of 
continuous incommensurate spiral orders 8 that yield a le- 
gal spectrum, but are unphysical at all. In contrast, the 
spectrum from LSWT has areas of negative energy for 
the collinear configuration when 1.2 < /x = J'/ J < 2.0, 
(X being the anisotropy defined below. But a variety of 
previous studies^— £ suggest that the collinear order is 
realized in this regime. It is believed that the instability 
of LSWT manifests the insufficiency of approximation. 
Therefore one of the objectives of this paper is to see 
if a NLSWT can fill this gap, with the negative energy 
excitation being considered in a novel way. 

Now let us clarify the model. J' mentioned above is the 
inter-chain spin exchange strength and J intra-chain. See 
FigJT] for the illustration. The nearest neighbors along 
the three orientations, 0°, 60° and 120°, are defined as 
(i,j)$- Si and S3 respectively correspond to the primitive 
translation vectors t\ and t 2 shown in Fig[TJ x direction 
is along t\ . The Hamiltonian is 

H = E J * E i s i s ! + s ? s j + s ? s !) w 

Here we consider the anisotropy J\ = J, J 2 = J3 = J'. 

The expected magnetic orders are signaled by the con- 
tour peaks at wave number Q's in the spin structure fac- 
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FIG. 1: Spin 1/2 TAFHM with the anisotropic nearest neigh- 
bor couplings J and J' (J' /J > 0) indicated in the figure. 
The primitive translation vectors are denoted by ri and T2. 
The periodic boundary conditions are assumed. 



needs infinite order to do so. In particular for S = 1/2, 
LSWT is regarded as an insufficient expansion for certain 
combination of spin configuration and coupling parame- 
ter. Our NLSWT goes beyond LSWT one order higher 
in 1/5*. It takes the following expansion 
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The Hamiltonian is rewritten in the K-space after the 
Fourier transformation as 



H — Hn + Hr> + Hr> + Ha 



(4) 
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(b) n = 1.3 
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where Hq is a constant contribution to the ground state 
energy (GSE) 
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—cosQ x + —^-cos—cos — - — (5) 



with L being the lattice size. On the other hand, H% is 
the quantum fluctuation bilinear with respect to magnon 
operators. It is 



FIG. 2: Contour plots of the pin structure factor by EPT for a Hi = ^ 

6x6 TAFHM (a) Isotropic when fj, — 1.0 and (b) anisotropic fe 
[i — 1.3. (a) refers to the 120° spiral spin configuration, (b) 
collinear. with 
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tor, Fig 121 from our numerical simulation 5 . Typically (a) 
Q = (4^,0), etc, for the f20° spiral configuration and 

(b) Q = ^0, -^J , etc, collinear. They are used to rotate 

the laboratory coordinate frame to unify the direction of 
the local spin in SWT. 

Since many work o 8 ' 10 ~ — have been done to formulate 
SWT for the spin 1/2 TAFHM, we simply start from the 
unitarily transformed Hamiltonian 

H ' = E J s E M + cos {q ■ ? ) (^J + sfsf) 

s (i,j) s 

+ sm (Q • 5) (SfSJ - SfS?)] (2) 

It is the Holstein-Primakoff transformation that trans- 
forms the above spin Hamiltonian to the bosonic operator 
(magnon) version. In practice it is approximated by ex- 
pansion in l/S. For example, LSWT is in the order of 

1/2 

0((f/S)~ 1/z )- The expansion restores the spin statistics 
in large S limit even in LSWT. However, for finite S, it 



Ak — — JcosQx — 2J cos—^-cos — + J cosk. 

T i ( Qx v^SQy , -, \ k x \/3ky 

+ J cos — cos — + i cos — cos (7) 

I 2 2 J 2 2 w 

cosQ x - I 
B k =- J cosk x 

T i ( Qx v3Qji \ k x \[%ky . . 

— J cos — cos — f cos — cos 8 

I 2 2 J 2 2 w 

For the moment we focus on Hq and H2, which are 
the only terms considered by LSWT. H3 and H4 are left 
for the later discussion. There is a handy transformation 
named after Bobolyubov to diagonalize H 2 . Its usual 
definition is as follows 

ak = u k b k + v k V_ k 

a k = U k°k + v kb-k (9) 
The coefficients, Uk and Vk, are the roots of the following 
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equations 



A\ — B\ is then expanded around k = Q as follows 



2 S(^)A fc 



5(A fc )gfc 



(10) 



Note that S (A K ) = Sign (A k ) is used in ([TO]). Since (fT0|) 
is in second order, one usually complements it with the 
following convention 



2 2 i 

w fc - v k = 1 



(11) 



to preserve the same commutation relationship of 
magnon. Then i?2, after the Bogolyubov transformation, 
reads as 



Ho 



(12) 



where 



4 - = J ■ 



{u-lf-l 



3/i 2 



* - — 
v V3 



from which it is seen that the stable region of the collinear 
states is 2.0 < u < oo in LSWT. (Both A k and B k are 
positive all over the if -space). LSWT starts to yield the 
negative energy excitation in the area surrounding those 
Goldstone modes when /i decreases below 2.0. It is obvi- 
ously unphysical, meaning that LSWT fails to account for 
the collinear spin configuration in this parameter regime. 
On the other hand, a variety of studies, either theoret- 
ical or numerical including our own work, suggest that 
the collinear configuration is realized when 1.2 < fj, < oo. 
Therefore it is desirable to investigate it in a way dif- 
ferent from LSWT. To this end, we here develop a NL- 
SWT which extends the Bogolyubov transformation to 
handle negative energy modes. We also introduce a self- 
consistent iteration scheme. 



e k = S(A k )S(Al-Bl) J\Al-Bl 



(13) 



defines the excitation spectrum. Since the ground state 
is the vacuum of Bogolyubov bosons, H 2 contributes to 
GSE as follows 



e 'gs = o E ( £fe _ Ak 



(14) 



Therefore GSE calculated within LSWT is 

Egs = E% s + E' GS (15) 



Now look at the property of e k . It can be easily shown 
that A k = B k at k = (0, 0), and A k = —B k at k = ±Q. 
Therefore they are always the zero-energy points, the so- 
called Goldstone modes. It is the characteristics of a 
system which spontaneously breaks the symmetry in the 
thermodynamic limit. In addition, the fisrt derivative of 
a linear dispersion at the Goldstone modes defines the 
magnon velocity, indicating that the excitation is indeed 
massless, another characteristics of the spontaneous sym- 
metry breaking. It can also be shown that there might 
be areas in K-space that have the negative energy exci- 
tation as /j, varies. For example, a collinear configuration 
when Q = (0, 27r/\/3) has 



A k — J (—1 + 2/i + cosk x ) 

I k x \/2tky 
B k = J \ ucos—cos — - — 



(16) 



II. A NOVEL NONLINEAR SPIN WAVE 
THEORY 

First of all, it is possible that the negative energy exci- 
tation exists even for SWT in infinite order if the starting 
spin configuration is unphysical. The emergence of neg- 
ative energy excitation could also be merely due to an 
insufficient approximation even if that spin configuration 
physically exists. The example here is the negative en- 
ergy excitation in LSWT for the collinear configuration 
when u < 2.0, which can be further divided into fi < uq 
and no < [i < 2.0 where uq is the physical PTP of the 
collinear configuration. The negative energy excitation 
might exist in the former regime for infinite order SWT, 
showing a sign of instability of such collinear spin order. 
While in the latter regime, one should be able to get an 
all positive spectrum in SWT with a sufficient approxi- 
mation, see below. 

Our idea of handling the negative energy modes is sim- 
ply to proceed like in the Dirac theory of positron 13 . 
We apply the particle-hole (anti-particle) transformation, 
b — > , with a boson normal-ordering, thereby formally 
filling up those negative energy modes and regard the 
filled states as a new physical vacuum. Since this for- 
mally involves an infinite negative energy, it is a sort of 
singular perturbation theory. When combined with mean 
field treatment of non-linear terms, we arrive at an iter- 
ative scheme of normal ordering and mean field approxi- 
mations to eliminate negative energy modes and improve 
on the physical vacuum which consists of zero-point mo- 
tion of all positive energy spectra. Or equivalently, we 
can make a singular choice, 



vl 



1 



(18) 



instead of (TTT|) . for negative energy modes k, which leads 
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to the particle-hole transformation, 



one needs to define the following mean-field parameters 



\Uu — V h 



(19) 



, and to a modification of ([12"]) as 



where C* = 1 for positive energy modes and —1 for neg- 
ative ones. 

In the following two sub-sections we first look at H4 
due to its relative simplicity and then #3. After the 
formulation, we analyze the collinear spin configuration 
first. It is a self-consistent iterative treatment of H4 
alone, since H3 vanishes for the collinear configuration. 
This is followed by the overall correction combining H3 
and H4 together for a general non-collinear spin config- 
uration and anisotropy. 



A. Quartic term 



The quartic term written in the magnon operator reads 



as 



where 



H 4 — Ha + Ha + Ha 



(21) 



h\ = y^ Js ro.s 



(Q ■ t\ ojoiojoj (22) 

1 — cos IQ-Sj. 
g (<'■"]"]". 



4- 11,11*11 , a , + a\a\ai(Jj + ajaiaiaj) (23) 

etc. After the Bogolyubov transformation with cither 
regular boson or new boson with a normal ordering, it 
can be written as 



H 4 = 6E 4 + 5H 2 + H 4 



(24) 



The first term is the Hartree-Fock correction to the 
ground state energy and the second the magnon self- 
energy. Since the last term describes two-particle scatter- 
ing processe, it only yields higher order 1/5 corrections 
compared to 8E4 and 6H2 . So we neglect it here. Before 
deriving the explicit form of 8E4 and SH 2 , we apply the 
Hartree-Fock decoupling to (|2"2")l . Due to Wick's theorem, 



in 



k 

= ( a l a j) = ^Z C *klkV, 



(25) 
(26) 



A = (o f Oj) = ^2 C^ k u k v k (27) 

k 

n=(a 2 l )=J2c* k UkVk (28) 
k 

To include the effect of anisotropy, one also needs to de- 
fine ms and As respectively as follows 



mi — y a l a jy = C k cosk x v 

1 k 



A 2 = {a iaj ) S2 = ^2c* k cos 

k \ 



u k v k (29) 



etc, where u k and Vk are defined either to give the reg- 
ular bosonic commutation or the new commutation for 
the Bogolyubov transformation depending on the sign of 
the excitation energy. One should notice that in the ex- 
tended stable parameter region all the excitation would 
be positive after the self-consistent iteration converges. 
Therefore the mixture of bosons and new bosons (they 
don't interact because they are defined for different wave 
numbers) will evolve during the iteration, in the con- 
verged state, only regular bosons survive. The mean-field 
parameters are functional of both Q and k. On the other 
hand, j k only reflects the average of the nearest neigh- 
bors in three orientations of a triangular lattice. It is 
defined as 



= - e lk ' s = - [ cosk x + 2cos—cos 
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(30) 



Now let us look at the following term which is a constant 
correction to the ground state energy, 



1 — cos 



(Q ■ S) (n 2 + m 2 + A 2 ) 
(Q • 6) 



1 + cos ( Q ■ S 



(2nA s + m s n) 



(2nm 5 + A 5 Q) 



(31) 



It emerges after the normal-ordering of the quartic term. 
The left-over two-operator terms are mean-field decou- 
pled, multiplied by the mean-field parameters. For ex- 
ample, 
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+ ^«l a ]) a i a j + ( a t a j) a \ a ) 



+ { a\a j ajdi + (a]ai ) ajcij (32) 



The correction to the quantum fluctuation of LSW from 
the quartic term reads in the K-sp&ce as follows 

SH 2 = SA k a\a k - ^5B k (a\<J_ k + a k a^ k \ (33) 



where 

SA k _ 
2 



1 - cosQ x 1 + cosQ x 
cosQ x m\ H si n \ cosk x 



- /' ( 2/gm' + -^^<> - (1 + f Q ) n ) f, 



1 — cosQ x 
2 

1 + cosQj 



Ai + m (1 - fq) A' 
mi - m(1 + /q) m' 



2 



(cosQj; + 2fif Q ) n 
cosQ x — 1 cosQ x + 1 



(34) 

fl — cosQ x Ai I cosk x 



(/ Q -l)n + i±^Sl-2/ Q A' 



COSQ^ - 1 COsQz + 1 

mi H Ai 



(35) 



with the auxiliary mean-filed parameters being defined 
as follows 



m'=J2C* k fkvt 

k 

A' = Y, C *kfkU k v k (36) 



t k = cos— cos — - — 
2 2 



B. Cubic term 

The cubic term i?3 in (jj) reads after the Primakoff 
transformation as 



H, 



E 



2SJ* 



(37) 



Under the Bogolyubov transformation, it becomes 
#3=^[ri (q,k-q,k) b\b\_ k b k 

kq 

+ T 2 (q, -k - g, k) b\bt 

where 

ijV2S\ 



q-k^k 



h.c. 



(38) 



Ti (1,2,3) 



T 2 (l,2,3) 



2! 



[Z\ («i + wi) (m 2 w 3 + V 2 V 3 ) 



+ Z 2 (u 2 + V 2 ) (U1M3 + VxVz) 

+ Z 3 (u 3 + v 3 ) {uiv 2 + viu 2 )} (39) 

iJV2S tr7 . w . 

= — - [Zi (Ml + Ul) (u 2 w 3 + v 2 u 3 ) 

+ Z 2 (u 2 + v 2 ) {u\v 3 + viu 3 ) 

+ Z 3 (u 3 + v 3 ) {uiv 2 + viu 2 )] (40) 

Above we wrote the expressions symmetric to the wave 
numbers. 2! in Ti is from the the permutation of 1 <-» 2 
in b\b 2 b 3 {b\b 2 b 3 ). And 3! in T 2 is from the permuta- 
tion among all three subscripts of b\b\b\ (bib 2 b 3 ). This 
symmetric form is somehow more suitable for the numer- 
ical integration over the first BZ, especially when there is 
cancellation of singular terms. Similar to the one fromii 
but (f38|) now includes the anisotropy as follows 



Z (Q, k, fi) =sinQ x sink x + 2a {h Q h k + g Q g k ) (41) 

. k x \^3ky 

h k =sin — cos 

2 2 



V3k 



v 



q k —cos — sin- 
■ y 2 2 

The effective Hamitonian reads as 

Heff = Hq + H 3 



(42) 



where Hq either refers to the LSW term H 2 or any unper- 
turbed Hamiltonian. If the anisotropy is set to vanish, 
(|39|) coincides with the ones i n 11 ! 14 . Following the stan- 
dard diagram perturbation method (See Appendix for 
details), the lowest order contributions are 



E 



|ri 



(k,u>) = -18^ 



uj — e q — e k - q + 
|r 2 (q,k)\ 



iO 

2 



iO 



to the normal self- energies, and 

,^T 2 (q,k)Tt(q,-k) 



E 



{k,u) = 6 
E>'^) 



^ uj + e q + e k+q - iO 



q 



T 2 (q,-k)T\ (q,k) 
- e k - q + iO 



(43) 
(44) 

(45) 
(46) 
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to the anomalous self-energies. However the anomalous 
term does not contribute to the spectrum in the same 
order O (S ) as the normal one's. The new poles of the 
magnon Green's function will then only reflect the quar- 
tic term and the normal self-energies from the cubic term. 
The real part of the poles determines the renormalized 
spectrum while the imaginary part determines the damp- 
ing rate. Explicitly, the new pole is as follows 



e k = e k + ixk 



4 



E 



2|ri (q,k)\' 



i8|r 2 ( 3 ,*or 



iO 



c k+q 



(47) 



Since Ti.2 are expressed in terms of u k and v k , it is de- 
sirable to write out T 2 2 explicitly in order to avoid the 
uncertainty of the sign of u k and v k . For example, 



III. SELF-CONSISTENT NSWT ANALYSIS OF 
COLLINEAR SPIN CONFIGURATION 



The purpose of this paper is to investigate the phase 
diagram of TAFHM by NLSWT, taking into account the 
negative energy excitation issue with higher order cor- 
rections. For the collinear spin configuration, we need to 
study the spectrum renormalization due to alone. In 
the last section, we extend the Bogolyubove transforma- 
tion to reformulate SH 2 in (|33j) . the contribution from 
H4 to the quantum fluctuation H 2 in LSWT. For the 
collinear spin configuration, the coefficients in (1331) arc 



5A k = (2 - 4/i) n + AuA' - 2m 1 + (2mi - 2n) cosk x 

+ n{2Q- Am') f k (51) 

SB k = (O - 2Ai) cosk x + n (4 A' - An) f k + A x - 2 urn' 

(52) 



with 



2rj(l,2,3) 
SJ 2 



2=1 



: (ui + vif (u 2 u 3 + v 2 v 3 ) 2 

5(Ai)(Ai + Bi) 
2V\Af-Bf\ 

( S (A 2 A 3 ) (A 2 A 3 + B 2 B 3 ) 
\^\{Al-Bl)(Al-Bl)\ 



(48) 



2°3 



(49) 



The correction from H4 to the ground state energy, (|31j) , 
is written explicitly for the collinear configuration as well 



6E 4 = (1 - 2a) (n 2 + r 
- (2nm + AO) 



2u (2nA + mfi) 



(53) 



Now we look at the spectrum renormalization due to 
5H 2 . One can obtain it by 



A k SA k -B k SB k 
e k = e k + S (A k ) : — : (54) 



or by 



and 

2T 4 
Z\Z 2 



--A ( Ul 



2°3 



f ui) (u 2 u 3 
S{A 1 ){A 1 - 



v 2 v 3 ) (u 2 
Bi) 



^J\A\-B\\ 
S(A 2 )(A 2 +B 2 



S(A 3 )(A 3 



B 2 2 \ 



B 3 



S(A 1 A 2 A 3 )(A 1 



S3 2 1 



v 2 ) (uiu 3 + VlV 3 ) 



B 1 ){A 2 + B 2 ){A 3 + B 3 ) 



Bj) 



(A 



BJ) 



(A 



B 2 3 ) 



(50) 



etc. T\ is formulated likewise. Note that the sign of the 
excitation energy of the unperturbed Hamiltonian is in- 
cluded. In the case where the unperturbed spectrum is 
positive, it has the same expression as after the usual Bo- 
golyubov transformation. However, (f50|) should be used 
whenever one considers NLSWT treatment involving a 
problematic spectrum from LSWT with a small area of 
negative energy. More discussion is given in the next two 
sections. 



€ k - Ct 



(A K + SA k ) 2 -(B K + SB k y 



(55) 



where is determined in the same way as . One im- 
mediately sees that when the correction 5 A k and SB k are 
small compared with A k and B k (l54|) is almost equiva- 
lent to ([S3)) except that (|M)) has singularity issue at some 
Goldstone modes where the sign of SB k / SA k is opposite 
to that of B k /A k . Physically this singularity will be can- 
celled either in the summation over the wave numbers 
in the first BZ or by the other source such as the cubic 
term in the same order as we will discuss later. However, 
it is challenging to numerically carry it out. Another 
drawback of (|54|) is that the mean-field quantities do not 
contain the information of the higher order correction. 
Therefore we prefer to use (|55p . In fact, the first way can 
be formulated from the Green's function using H 2 as the 
unperturbed Hamiltonian. The second way can be ob- 
tained by absorbing 5H 2 into the unperturbed Hamilto- 
nian. We believe that when the correction is small, these 
two approximations shall both yield satisfactory results. 
However so far they are only a half-cooking recipe, be- 
cause even if they can obtain positive spectrum, neither 
guarantees the non-recurrance of negative energy prob- 
lem. The situation would be more sound if a convergence 
can be provided by a self-consistent iteration. Below we 
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FIG. 3: LSWT and NLSWT results for the collinear spin 
configuration on a square lattice by approaching fi to oo and 

setting Q= (0,||). 



(a) /i = 1.2 (b) fi = 1.3 
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FIG. 4: Comparison between the contour plots of the exci- 
tation spectrum in the first BZ for (a) /j, = 1.2 and (b) 1.3. 
Both spectra are calculated by NLSWT with a converged self- 
consistent iteration. The dark blue region in (a) refers to the 
negative energy excitation, (b) has no negative energy exci- 
tation. 



illustrate the idea. 

First we start the self-consistent iteration from the 
ground state of LSWT to get the mean-field parame- 
ters. If it has the negative energy excitation for some 
wave numbers, it is not the lowest energy state but a 
meta-stable state. When we calculate the mean-field pa- 
rameters, the new commutation ()18j) should be used to 
collect the contribution of those wave numbers associated 
with the negative energy excitation. At the end of each 
iteration one uses the obtained mean-field parameter to 
update the excitation spectrum. In other words, one re- 
assesses which excitation is regular boson or the new bo- 
son in the first BZ. In the next iteration, the ground 
state is chosen to be the vacuum of those updated bo- 
son and new boson. Still it is a meta-stable state if the 
negative excitation does not vanish in the whole spec- 
trum. This iteration is carried on until convergence. In 
the end, there are two possibilities. First, the excitation 
spectrum becomes all positive, and the ground state has 
the lowest energy, a physical ground state. Second, iter- 
ation is converged but there still remain some negative 
energy excitations, signaling the instability of assumed 
magnetic order. 

We perform our calculation on a 1800 x 1800 lattice for 
which the finite size effect can be neglected. We want to 
first test the validity of our formulation. To this end, we 

set Q = (o, and fi — > oo to check the well known re- 
sult of a square lattice with a collinear spin configuration. 
In this case, there is no negative energy excitation issue 
even in LSWT. Indeed, in FigEJone of our NLSWT calcu- 
lations, which uses (fMf and doesn't take a self-consistent 
iteration process, agrees well with-^. We also see in 
Fi^3]that the energy after the self-consitent iteration con- 
verged is obviously further from the believed ground state 
energy shown by star due to variational Monte Carlo^— 



that is tested by our own nemerical simulation^. Since 
self-consistency means to take higher order correlations 
into account, NLSWT with self-consistent iteration is ex- 
pected to give a better ground state energy than without 
self-consistent iteration. That it is not the case reflects a 
fact that SWT is not variational and hence an approach 
to the correct value need not be monotonic. 

On the other hand, the self-consistent iteration works 
better in the intermediate parameter regime in the sense 
that it improves on the excitation spectrum in 1.3 < 
/j < 2.00 and that it avoids an artificial oscillation of 
the NLSWT results without a self-consistent iteration. 
Furthermore, the discontinuity at fi « 1.3 is consistent 
with the fact that the improvement of spectrum takes 
effect when /i > 1.3 as shown in FigJU That is to say, 
it provides a clear and consistent sign that there is a 
phase transition from the collinear spin configuration to 
another configuration when fj, decreases below 1.3. 

It is known that PTP of collinear configuration in 
LSWT is [i — 2.0. In contrast, it is /i w 1.3 in our NL- 
SWT. So we compare the excitation spectrum between 
LSWT and our NLSWT around these two points. Figj5] 
shows the comparison (a) around the PTP of NLSWT 
and (b) around the PTP of LSWT. Besides the obvi- 
ous occurance of negative energy excitation when fi de- 
creases below their PTPs, the softening of the spin wave 
at T-point deserves an emphasis. The strongest softening 
happens when the sound velocity is calculated along the 
direction from T-point to X-point. For example, it is zero 
at /j, = 2.0 in LSWT, implying a breakdown. 

It is interesting to see how PTP evolves during the iter- 
ation. We change fi in the range of 1.2 < [i < 2.0 and for 
each value investigate the excitation spectrum at every 
step of iteration. We are then able to determine PTP 
vs iteration steps (IS) shown in Fig|6l One sees that 
the PTP of the collinear spin configuration approaches 
quickly to the convergence right after a few iteration. It 
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(a) Around NLSWT PTP (b) Around LSWT PTP 



5 



FIG. 5: Comparison between the spectra of LSWT and NL- 
SWT around the phase transition point (PTP) at (a) fi = 1.3 
for NLSWT and (b) 2.0 for LSWT. 
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FIG. 7: The finite-size scaling of the contribution to the 
ground state energy from (a) CT and (b) QT. This scaling 
is for the 120° spiral configuration on an isotropic lattice 
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FIG. 6: The star (blue online) show PTP of LSWT. It moves 
to n = 1.52 immediately after one iteration. Eventually it 
converges to /i « 1.3. 



is an important feature, not only because it reveals that 
the new commutation relation renormalizes the spectrum 
systematically (In fact, without a correct treatment of the 
negative energy excitation we found no way to renormal- 
ize the spectrum to the physically meaningful one while 
the correct symmetry, i.e., the Goldstone mode, is in- 
tact), but also because this nice tendency yields a useful 
information about the phase transition even when it does 
not converge. 



IV. NLSWT ANALYSIS OF NON-COLLINEAR 
SPIN CONFIGURATION 

We also want to test our method for non-collinear spin 
configurations and anisotropy by applying it to the limit- 
ing case which can be readily compared with the existing 
calculations. We discuss in Sec lIV A] the commensurate 



spiral order, Q = (47r/3, 0) and fi= 1. It is a 120° three 
sub-lattice spin configuration on the isotropic lattice. In 
Sec lIV Bl we extend the discussion to the incommensurate 
spiral order. 



A. 120° Commensurate Spiral Order 

Note that the contributions from the cubic term (CT) 
and from the quartic term (QT) both have the singularity 

at k = ±Q. From Figj7]one sees that their singularities 
exactly cancel out only at the infinite limit. Since the 
numerical calculation can be done only on a finite lat- 
tice, a finite-size scaling is necessary to extrapolate the 
precise corrections. For example we show two quantities 
that are related to the contribution to GSE, 13 from CT 
in (a) of Fig(7] and 14 from QT in (b) of the same fig- 
ure. Particularly for an isotropic lattice 14 = —88E4/3J 
with 5E 4 being defined in JS) and 13 = —45E 3 /3J with 
5E3 being the contribution to GSE from CT. 13 reads as 
follows 



13 = 



, |r 2 (g,fc)| 2 



(56) 



One sees from FigJT] that 13 converges much slower than 
14. Note that our direct numerial simulation is much 
heavier compared to the finite-size extrapolation using 
clusters with different aspect ratios 1 ^, and gives a clear 
convergence as seen in FiglT] The extrapolated values of 
13 and 14 agree wit h 11 ' 14 denoted by blue stars. In addi- 
tion, it is well known that in a direct numerical simulation 
a small number S, comparable with the increment of nu- 
merical integration in the K space, should be inserted in 
the denominator of any fractional number. It is to pre- 
vent the overflow in the numerical simulation. Therefore 
finite-size scaling using different S is presented in the fig- 
ure. Smaller 6 gives a faster convergence. The standard 
choice 5 — tt/L is also good as seen in the figure. 
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FIG. 8: The finite-size scaling of the spectrum for the 120° 
spiral configuration on an isotropic lattice. The slice from T 
point, k = (0, 0) to K point, k = Q — (47r/3, 0), is shown. 



We then use 8 — n/L to perform the numerical simu- 
lation for the spectrum. A finite-size scaling is shown in 
FigJS] It is clear that the singularities from CT and QT 
at if-point indeed exactly cancel out at the thermody- 
namic limit but it scales very slowly. The spectra from 
top to bottom in the figure account for different lattice 
sizes, from L = 300 to 76800. The extrapolated spectrum 
agrees withii. 

It is also notable that there exists discontinuity in the 
spectrum, shown as X-point in Fig|7] It could be cured 
by some techniques. For example the damping can be 
absorbed into the Dyson equation to further renormalize 
the spectrum like in the off-shell approximation 11 . How- 
ever we don't go in this direction for two reasons. First, 
the 120° spiral order when \i = 1.0 does not have nega- 
tive energy excitation despite the discontinuity. Second, 
those further treatments like off-shell approximation are 
unlikely to cure the negative energy excitation when fi 
varies away from 1.0 while the configuration is fixed. In- 
stead, in the next section we discuss the instability of the 
incommensurate spiral order in NLSWT, which is origi- 
nally stable in LSWT at specific values of /i. 



B. Incommensurate Sprial Order 
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FIG. 9: The finite-size scaling of the spectrum for incommen- 
surate spiral configurations. A slice is shown, from V point to 
if -point of (a) (7tt/6, 0) and (b) (tt/2, 2tt/v / 3) 



into < fi < 0.27. The difference can be explained if the 
reduction, the difference between the magnetization and 
the classical value (1/2 for spin one half), is considered. 
It is found that the reduction exceeds the classical value 
when < /j, < 0.27 in LSWT, indicating its breakdown 
in this regime and that the state is disordered. 

We study these incommensurate spiral configurations 
by NLSWT in this session. A few wave numbers are 
chosen to carry out the numerical simulation. The cor- 
responding value of fx resides on both sides of 1.0. And 
only two typical spectra are plotted. They are (a) of 
FigHwhen Q = (7tt/6,0) and (b) (tt/2, 27T/V3) . The 
spectrum is plotted along the same slice as in FigJHl with 
only the vicinity of X-point, one of the Goldstone modes, 
being shown. 

The solid lines show the extrapolated spectrum by the 
finite-size scaling. They exhibit both the negative energy 
excitation and a discontinuity at if-point. The relation- 
ship between the distance from the commensurate wave 
number, AQ = Q x — 4ir/3 with Q x being the unfolded 
incommensurate wave number, and A, the depth of the 
negative energy, i.e., half of the kink since the disconti- 
nuity is symmetric about 0, is plotted in Fig[TU] It is 
linear near the commensurate spiral wave number. Al- 
though only a few wave numbers are checked, it seems 
in Fig llOl that all the incommensurate wave numbers are 
corresponding to spectra with negative energy excitation 
and discontinuity at if-point. Namely, the incommensu- 
rate spiral spin configuration is unstable in the NLSWT. 



There exist incommensurate spiral states in LSWT 
when /i varies away from 1.0. The relationship be- 
tween /i and the ordering wave number is: Q = 
(2COS- 1 (-m/2) ,0) when 0.27 < (i < 1.0 and Q = 
(2ir - 2ms- 1 (-/i/2) , 2tt/V3) when 1.0 < fi < 2.0. Note 
that the splitting into branches is just due to the folding 
into the first BZ of an originally continuous line. A sim- 
ple calculation shows that this relationship is identical 
to the one in the classical picture except that it extends 



C. Summary and Proposal of a Self-Consistent 
Iteration 

From the last two sections we conclude that in this 
NLSWT only the commensurate 120° spiral order is sta- 
ble among all the spiral configurations. But it would be 
interesting to see if special treatments like the off-shell 
approximation can cure the discontinuity at the if-point, 
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FIG. 10: The finite-size scaling of the spectrum for a 120° 
spiral configuration on an isotropic lattice. The slice from T 
point, k = (0, 0) to K point, k = Q — (47r/3, 0), is shown. 



say, at least for some incommensurate spiral wave num- 
ber. If so, the depth of the negative energy will then 
disappear because it is tightly associated with the discon- 
tinuity However, this perspective is essentially different 
from what off-shell approximation does for the commen- 
surate spiral state where the discontinuity at A-point is 
smoothed by it while the spectrum at A'-point is not al- 
tered. 



V. CONCLUSION AND DISCUSSION 

We have extended the Bogolyubov transformation by 
considering a new commutation relation for the bosonic 
quasi-particle which has negative energy. This new trans- 
formation, equivalent to a boson normal ordering, bridges 
a self-consistent iteration from a spectrum with areas of 
negative energy excitation in LSWT to an expected posi- 
tive energy spectrum in NLSWT. The new stable regime 
1.3 < (J, < oo for the collinear spin configuration on a 
spin 1/2 TAFHM is comparable with the ones by other 
methods 4 ^— Considering the clear physical picture of 
both the ground state and the excitation spectrum pro- 
vided by SWT, such improvement of the approximation 
precision over the previous SWT studies ^ 10 ' 12 stimulates 
our interest to apply the same scheme to non-collinear 
spin configurations. Indeed our study shows that among 
spiral configurations, only the 120° commensurate sprial 
is stable at fi — 1.0 while all incommensurate ones are un- 
stable. This findling has two-fold meanings. First, most 
incommensurate sprial states are not seen by previous 
studies, which coincides with this finding. Second, how- 
ever, certain incommensurate spiral state is seen. For in- 



stance, Q = (n/2, 27r/-\/3), namely a four sub-lattice spi- 
ral state has been seen in a very narrow window around 
fi = 1.18 by our recent Eangtanglement Perturbation 
Theory study*. It implies that due to the approximate 
nature of SWT, especially when the spin is as small as 
1/2, NLSWT still misses certain details in the phase dia- 
gram. Moreover, our more elaborate self-consistent itera- 
tion study for a non-collinear spin configuration does not 
give a satisfactory result of an expected broadening of the 
stable region for 120° spiral state. A further development 
of NLSWT along this line shall be highly desirable. 
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Appendix A: Self-Energy 

In this appendix, we show the formulation of the self- 
energy in the lowest order for CT, for instance, ^2 (k, w). 
The green function is defined as follows 

G (k, r) = - (T T b k (r) b{ (0)) = - (T T U (0) b k (r) b{ (0)) 

(Al) 

with 

UW) = J2(-r) ~ dT 1 ---dr n (T T H'(T 1 )---H'(T n )) 

(A2) 

The perturbative Hamiltonian is H' = H 4 + H 3 . One 
immediately sees that iJ 4 , if one choose not to absorb it 
to H°, contributes in (|A2[) as the first order, it will give 
e k , just a HF approximated term as we discussed before. 
Now we focus on the self-energy due to H 3 , which arises 
from the expansion of (|A2|) in the second order: 

c/(2) (/3) = h h. jf dr ^ H " ^ h3 ^ ( A3 ) 

Therefore 

G (2) (k, r)=-~ jf" dndr 2 (H 3 ( Tl ) H 3 (r 2 ) b k (r) b{ (0)) 

(A4) 

Using ([55)1 we have 
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H 3 ( n ) H 3 (r 2 ) = r i *l " <&> fc i) E r * fc2 " ^ fc2 ) 6 * ( r i) ( r i) 6 *> ( T i) 6 * ( T2 ) 6 »-te ( T2 ) b L ( T2 ) + ' 

gifci 92^2 

(A5) 



The first term above will give rise to ^ (k,oj), Now we orem, one of terms in the form of operator-pair to be 
show the formulation step by step. Due to Wick's the- summed in (|A4p is 



(Trbl (n) 6, 2 (r 2 )) (iV^ ( Tl ) 6 fc2 _ 92 (75)) (T T b kl (n) b{ (0)) (t t &[, 2 (r 2 ) 6 fc (r)) 

= (T T b q2 (r 2 ) &J t (n)) 5 91)52 (T T b k . 2 - q2 (r 2 ) (n)) 4 1 - gi ,fc 2 - g2 (^-r&fc! Oi) &[. (0)) 5 fcl , fc {T T b k (r) b\ 2 (r 2 )) 4,fc 2 

=G (<7i,r 2 - Ti)(5, Zl! g 2 Go (fci - #1,7-2 - Ti ) 4 1 - 9l ,fc 2 -g 2 G (fc,ri) 5fc 1; fcG (fc,r - r 2 ) 6 fcjfc2 (A6) 

I 



Where Go is the unperturbed Green function, 
G (fc,r) = -(T T b k (t) fct (0)) 



(A7) 



Substituting (|A7[) into (| A6|) and using the kronig-delta 
function one has 



Gi 2) (k,r) 



1 

2ft 2 



dridr 2 E l r i (?) k )\ 2 G o (l> T 2 - n) 



Go (k -q,T 2 - n) G (k,r- r 2 ) G (A,n) 

(A8) 

The first subscript stands for the lowest order while the 
second for only part of terms. After the Frourier trans- 
formation, (|A8p reads as 

G< 2) (k,iu) 

=4 E E i r i («. fc )i 2 G " («> G ° ( fc ' ^) 

q u' 

Go (fc, Go (fc — g, iu> — ioj') 
=G n {k,iu)Y^ 1 (k,uj)G (k,iuj) (A9) 

Thus one gets the implicit expression of part of the first 
normal self-energy 

1 



E^^EE^^r 



or q 



G {q,iuj')G {k-q,iuj-iLo r ) (A10) 

If one needs its explicit expression, the unperturbed 
Green function above must be integrated out with re- 



spect to q and u>'. Take the definition^ 
G (k.iuj) =- — 

ILU — U! k 

And note that the boson distribution function is 

1 



n B (z) 



a0Z _ 1 



(All) 



(A12) 



Its poles are the Matsubara frequencies for boson. One 
then transforms the frequency summation to an integral 



B*-> '2m L 
Explicitly, one defines a contour integral as follows 
I 

=^E l Fl ^,fc)| 2 j^ dzn B (z 



(A13) 



1 



1 



=0 



(IA14I) can be rewritten as follows 

I = R 1 +R 2 =0 



Z — IU1 + Cd k - q Z — UJ q 

(A14) 
(A15) 



where R\ — — Yll {k,w), the residue on poles of ns, 
and i? 2 is the residue on poles of Go's. After analytic 
continuation, the expression of (k, w) is 
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£° (*, w) = S (g 9 ' fc)|2 ( «b -— " =7 + [1 - n B {oJk- q )] ■ r 1 } (A16) 

^1 ^ 2 [ W g + CJfc_ g - UJ - lS Ul + Id - 0J k - q - LU q J 



r 



At the zero temperature, tib — 0, its expression is as fol- 
lows 

V a (fc^) = lV |ri(g ' fc)|2 - (A17) 
' 2 ^ LJ-e q -e k - q +iS v ; 

The other terms due to Wick's theorem can be obtained 
similarly. The final expression of J] a (k,co) is 

r(^) = 2V |ri(g ' fc)|2 - (A18) 
^ ^u-e q -e k - q + i6 
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